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16 0+ 16 nature of the superdeformed band of 32 S 
and the evolution of the molecular structure 
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The relation between the superdeformed band of 32 S and 16 0+ 16 molecular bands is studied 
by the deformed-base antisymmetrized molecular dynamics with the Gogny D1S force. It is found 
that the obtained superdeformed band members of 32 S have considerable amount of the 16 0+ lfi O 
component. Above the superdeformed band, we have obtained two excited rotational bands which 
have more prominent character of the 16 0+ 16 molecular band. These three rotational bands are 
regarded as a series of 16 0+ 16 molecular bands which were predicted by using the unique ie O - 
16 O optical potentil. As the excitation energy and principal quantum number of the relative motion 
increase, the 16 0+ 16 cluster structure becomes more prominent but at the same time, the band 
members are fragmented into several states. 

PACS numbers: Valid PACS appear here 



The properties of the 16 0+ 16 molecular bands have 
been studied by many authors with the microscopic clus- 
ter models for many years Q. Despite of these studies, 
the microscopic models have not been able to give a con- 
clusive answer. One of the reasons is the fact that the 
number of the molecuar band, the excitation energy of 
the band head, and the moment of the inertia strongly 
depend on the effective nuclear force. Recently, a rather 
conclusive answer was given by the studies with the 
macroscopic model 0, U In those studies, the authors 
used the unique potical potential for the 16 0- 16 system 
which was determind in 1990's |4| after the first discov- 
ery of the nuclear rainbow in 1989 [5j- These studies gave 
the following answers for the lowest three rotatinal bands 
whose principal quantum numbers N=2n+L of the rela- 
tive motion between clusters are N=24, 26 and 28, respec- 
tively: The lowest Pauli-allowed rotational band (N=24) 
starts from the + state located at about 9 MeV in 
the excitation energy (about 8 MeV below the 16 0+ 16 
thrcthold), and the energy gap between the N=24 and 
N=26 bands and that between N=26 and N=28 bands 
are both approximately 10 MeV. In Ref . 2] , it is proposed 
that the observed 16 0+ 16 molecular states correspond 
to the third band whose principal quantum number is 
N=28. 

Besides the cluster models, in these days, the superde- 
formed structure of 32 S has been studied by many authors 
with the mean- field theories 0, 0, • ^ is largely because 
the superdeformed structure of 32 S is regarded as a key to 
understand the relation between the superdeformed state 
and the molecular structure. Indeed, by the HF(B) cal- 
culations, it is shown that the superdeformed minimum 
of the energy surface is well established in each angular 
momentum, and at the superdeformed local minimum, 
the wave function shows the two-center-like character. It 
is also notable that many of the mean-filed calculations 



predict that the superdeformed band starts from the + 
located at around 10 MeV which agrees with the band 
head energy of the N=24 band obtained from the unique 
optical potential. Therefore, it is enough conceivable that 
the superdeformed band obtained by the mean-field cal- 
culations and the lowest Pauli-allowed 16 0+ 16 molecu- 
lar band (N=24) are identical. 

In the present study, we aim at clarifying the rela- 
tion between the superdeformed state and the 16 0+ 16 
molecular structure. The objectives of the present study 
are summarized as follows. (1) To what extent the su- 
perdeformed state and the 16 0+ 16 molecular structure 
are related? : In the unique optical potential analysis, 
the factors which distort the 16 0+ 16 cluster structure 
such as the effects of the spin-orbit force and the for- 
mation of the deformed one-body field are not treated 
directly. Instead, these factors are renormalized into the 
optical potential through the extrapolation to the low- 
energy region. When one treats these factors veraciously 
by the microscopic models, the pure 16 0+ 16 cluster 
structure will be distorted and will have a deformed one- 
body field character. In other words, the superdeformed 
states in the mean-field models and the states of the low- 
est Pauli-allowed 16 0+ 16 band of the unique optical 
potential will be the states which have both characters 
of the deformed one-body filed structure and two cluster 
structure. Therefore, it will be important to study to 
what extent each state has the 16 0+ 16 character and 
the deformed one-body field character. (2) Do the excited 
states exist in which the excitation energy is spent to ex- 
cite the relative motion between the clusters? Do they 
correspond to the N=26 and 28 bands obtained from the 
unique optical potential? : When we believe that the su- 
perdeformed states of 32 S have the considerable 16 0+ 16 
component, we can expect the excitation mode in which 
the excitation energy is used to excite the relative motion 
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between the clusters. 

The deformed-base antisymmetrized molecular dy- 
namics (deformed-base AMD) [l| combined with the gen- 
erator coordinate method (GCM) has been used with the 
Gogny D1S force [Io|. For the sake of the selfcontained- 
ness, we briefly explain this framework. For more details, 
the reader is refered to the references ^]J. The intrinsic 
basis wave function of the system is expressed by 
a Slater determinant of single-particle wave packets tpi. 
Each single-particle wave packet is composed of spatial 
part 4>i, spin part Xi an d isospin part r,;. The spatial part 
has the form of the deformed Gaussian centered at Z^. 

®int = ^=det{^(ij)}, (1) 

Vi(rj) = M r j)Xin, (2) 
4>i{r) = cxp{-^(x - Z„) 2 - u y {y - Z iy ) 2 

-v z (z-Z. lz ) 2 }, (3) 
X« = a iX] + PtXii Ti = proton or neutron. (4) 

Here, the centroids of the Gaussian Zj and the spin di- 
rection on and Pi are complex parameters and are de- 
pendent on each particle. The width parameters v Xl v y 
and v z are common to all particles. These variational 
parameters {Z i: a i; /3j and (v x ,v y , v z )) are determind 
by the variational calculation. The variational calcula- 
ton is made after the parity projection by using parity- 
projected wave funciton ^ = ^mt as the varia- 
tional wave function. In this study, the variational cal- 
culation is made under the constraint on the nuclear de- 
formation parameter (3. The advantage of the deformed 
Gaussian basis as the single-particle wave packet is that 
it is possible to describe both the deformed one-body- 
field structure and the cluster structure as well as their 
mixed structure within the same framework. We can con- 
firm this feature when we consider the two limits of the 
nuclear structure described by this wave function, the 
deformed- harmonic-oscillator limit and the cluster limit. 
The deformed-harmonic-oscillator limit is reached when 
the centroids of all single-particle wave packets (Re Zj) 
are at the center of the nucleus and the single-particle 
wave packets are deformed. On the contrary, the cluster 
limit is obtained when the centroids of the single-particle 
wave packets are separeted into the centers of the con- 
stituent clusters and the single-particle wave packets are 
spherical. 

After the constrained variational calculation for 
we superpose the optimized wave functions for various 
deformation parameters (GCM): 

^ = cP J MK ^{(3 ) + c'P J MK ,^{(3' ) + (5) 

where $ 7r (/3o) is the optimized wave function under the 
constraint of the nuclear deformation parameter = f3o 
and P'l IK is the angular momentum projector. The 
cofhcients c, c',... are determined by the diagonalization 




nuclear deformation p" J 

FIG. 1: The energy surface as a function of the nuclear de- 
formation j3 for 32 S. Dashed curve is for the J = state 
obtained by the 16 0+ 16 Brink wave function. Solid curves 
are for J=0, 2, 4, 6 and 8 states obtained by the deformed-base 
AMD. The matter density distribution of the deformed-base 
AMD wave function at the superdeformed minimum (J=0) is 
also shown. Small black circles in the density distribution rep- 
resent the centroids of the single-particle wave packets ReZi 

of the Hamiltonian. When the intrinsic wave function 
$ 7r (/3 ) has a prominent cluster structure, the superpo- 
sition of the wave functions will improve mainly the de- 
scription of the relative wave function between clusters, 
while it will improve mainly the description of the mean- 
field and other correlations when $> v (f3o) has a prominent 
mean-field character. 

Below we explain how we estimate the amount of the 
16 0+ 16 component in the obtained wave function <& j7r . 
We rewrite the $ j7r by decomposing it into the 16 0+ 16 
component and the residual part, 

^ = a^{ X j(r)Y JO (f)0( 16 O)<K 16 O)} + V~l - a 2 $„(6) 

where A is the antisymmetrizer, r is the relative coordi- 
nate between two 16 clusters and cf>( 16 0) is the internal 
wave function of the 16 cluster. x,/(r) is so normal- 
ized that A{x.j{r)Yjo{f)(l)( w O)(j>( i6 0)} is normalized to 
unity. $ r is orthogonal to the 16 0+ 16 cluster space. 
It is possible to evaluate the amount of the 16 0+ 16 
component w J = \a\ 2 = |($ J,r |P/|$ j7r )| 2 with the pro- 
jection operator P[ = J2 a l<?a)(<?al onto the 16 0+ 16 
cluster model space, where <p^ is the orthonormalized set 
of the spin projected Brink wave functions for 16 0+ 16 
system. In this study, 26 spin-projected Brink wave func- 
tions with the inter-cluster distance ranging from 0.5 fm 
to 13.0 fm are orthonormalized to construct <p J a . 

Before discussing the results of the deformed-base 
AMD+GCM, we show the result of the calculation in 
which we use the intrinsic wave functions which have the 
pure 16 0+ 16 configuration. The Brink wave function 
is used for the intrinsic 16 0+ 16 wave function. There- 
fore, this result is the same as the traditional cluster 
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model calculation of 16 0+ 16 RGM (resonating group 
method) and 16 0+ 16 GCM. In Fig. |TJ the energy sur- 
face for the J = state as a function of the nuclear 
deformation is shown (dashed curve). We note that the 
16 0+ 16 configuration describes more than Atuv excited 
states relative to the ground state of 32 S, therefore the 
ground state and the normal deformed states are not in- 
cluded in this energy surface. It has the energy minimum 
at /3 = 0.73 (inter-cluster distance is 5.0 fm) which cor- 
responds to two touching 16 O. The minimum energy is 
about 8 MeV higher than the 16 0+ 16 threthold en- 
ergy. After the GCM calculation along this energy sur- 
face, we have obtained three rotational bands which have 
the quantum numbers of the relative motion N=24, 26 
and 28, respectively, and are shown by dashed lines in 
FIG. [21 However, their energies are too high to coin- 
cide with the rotational bands obtained from the unique 
optical potential and also with the superdeformed band 
obtained from the HF(B) calculations. The energy gaps 
between these bands (4 MeV between N=24 and N=26, 
and 6 MeV between N=26 and N=28 in the case of + 
states) are much smaller than the results from the unique 
optical potential. We think that these deviations come 
from the fact that the effects which distort the cluster 
structure are neglected in the Brink wave function. We 
will see below that in fact the effects of the distortion are 
fairly large. 

In Fig. ^ the energy surfaces for J = to J = 8 states 
obtained by the deformed-base AMD+GCM calculation 
are also shown. Since the AMD wave function does not 
assume any cluster configuration, the normal deformed 
states also appear in these energy surfaces. In the nor- 
mal deformed states, the prolately deformed states and 
the triaxially deformed states (7 = 6° ~ 30°) are en- 
ergetically degenerate. After the GCM calculation, pro- 
late wave functions mainly contribute to the ground band 
while the triaxially deformed wave functions to the first 
excited band. Their excitation energies and the intra- 
band E2 trainsition probabilities (Table HJ) show reason- 
able agreement with experiments and are consistent with 
the results of the HFB calculation with the Gogny D1S 
force |3] , though the total binding energy of the ground 
state underestimates the experimental data by about 2.0 
MeV. 

We also find that the behavior of the energy surface 
around the excited local minimum is similar to that of the 
HF(B) calculations. In each angular momentum state, 
the superdeformed mimimum is well developed and the 
+ excitation energy relative to the normal deformed 
ground state is around 10 MeV. The energy difference be- 
tween the deformed-base AMD and the 16 0+ 16 calcu- 
lation at the superdeformed minimum is about 10 MeV, 
which implies a fairly large effect of the distortion of 
the cluster structure. Indeed, the deformed-base AMD 
wave function deviate from the pure cluster limit. At 
the superdeformed limit, the single-particle wave pack- 
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FIG. 2: The excitation energies of the N=24, 26 and 28 
band members obtained by the 16 0+ 16 Brink wave func- 
tion (dashed lines) and the deformed-base AMD+GCM (solid 
lines). The N=26 and N=28 band members are fragmented 
into several states in the deformed-base AMD+GCM calcula- 
tion and the averaged energies Eav are shown for these bands. 



ets are prolately deformed (v x = v y = 0.160 fm - and 
v z = 0.115 fm -2 ), and the distance between the centroids 
of the single-particle wave packets are rather small (3.1 
fm), though they are separated into two parts (centers of 
two 16 0) and show two-center nature. The energy gain 
by the deformed-base AMD function compared to the 
pure 16 0+ 16 wave function mainly comes from the two- 
body spin-orbit force and the density dependent force. 
In the case of the deformed-base AMD wave function, 
the expectation value of the two-body spin-orbit force is 
about -4.5 MeV which must be zero in the pure 16 0+ 16 
wave function. Though its value is not so large, it lowers 
the excitation energy of the superdeformed state from 
the 16 0+ 16 limit. Compared to the 16 0+ 16 wave 
function, the expectation value of the repulsive density 
dependent force is about 6 MeV smaller in the deformed- 
base AMD wave function and it also indicates the non- 
small deviation from the 16 0+ 16 structure. Though 
the kinetic energy does not much contribute to lower the 
energy of the superdeformed band, its nature is different 
from the case of the pure 16 0+ 16 structure. At the 
superdeformed minimum, the single-particle wave pack- 
ets are prolately deformed and since the kinetic energy 
is almost linear to width parameter i>, the kinetic nergy 
to the z direction is eased. However, we have found that 
the wave function of the deformed-base AMD (which is 
parity and the angular momentum projected) at the su- 
perdeformed minimum still has a considerable amount of 
the 16 0+ 16 component, w J=0 = 0.57. 

After the GCM calculation by superposing the 
deformed-base AMD wave functions along the energy sur- 
face, we have obtained three rotational bands. The low- 
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TABLE I: Excitation energy Ex [MeV] and the intra-band E2 transition probabilities B(E2;J — » J — 2) [e 2 fm 4 ] of the ground 
band and the first excited band. 





ground (Theor) 


gr 


ound (Exp) 


band I (Theor) 


band I (Exp) 
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Ex 


B(E2) 
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B(E2) 


Ex 


B(E2) 


Ex 


B(E2) 













3.9 




3.778 
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2.3 


66 


2.23 


60±6 
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31 


4.282 
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5.75 


109 


4.459 


72±12 


10.1 


88 


6.852 


35.4^ 


6 


10.2 


130 


8.346 


> 22.2 


12.9 


98 


9.783 



deformed base AMD+GCM 
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FIG. 3: The excitation energies of the N=24, 26 and 28 band 
members obtained by the deformed-base AMD+GCM (solid 
lines) and the deformed-base AMD+( 16 0+ ls O)+GCM (dot- 
ted lines). The N=26 and N=28 band members are frag- 
mented into several states in both calculations and the aver- 
aged energies Eav are shown for these bands. The deforemd- 
base AMD+GCM results in this figure are the same as those 
of FIG.H 



est superdeformed band, the second lowest band and the 
third band have large 16 0+ 16 components with N=24, 
26 and 28, respectively. Thus these three bands corre- 
spond to the N=24, 26 and 28 bands of the unique opti- 
cal potential. We refer to these bands simply as N=24, 
26 and 28 bands (Fig. |2J). The band members of the 
N=26 and 28 bands are fragmented into several states. 
Namely, there are several states which have the 16 0+ 16 
component of the same principal quantum number N. 
For example, the + state of the N=26 band is frag- 
mented into three states in which w J=0 are 0.34, 0.13, 
and 0.25, respectively. In the following, we only show 
the averaged energies of these fragmented states to in- 
vestigate the gross feature of the band structure and the 
fragmentation is discussed later. The averaged energy is 
calculated by multiplying the overlap to* as the weight; 
El v = {w{E{ + w^E^ + ...)/{w{ + wl + ...). The exci- 
tation energies of the N=24 band members are lowered 
by about 3 MeV compared to those of the band formed 
by the states at the superdeformed energy minimums. 



This may be due to the improvement of the description 
of the mean-field and the collective motion by the GCM 
calculation, since in the N=24 band, the amount of the 
16 0+ 16 component is decreased by the GCM calcula- 
tion; w J=0 — 0.57 for the deformed-base AMD wave func- 
tion and w J=0 = 0.42 for the deformed-base AMD+GCM 
wave function. It is important that the energy gain of the 
deformed-base AMD+GCM band compared to the pure 
16 0+ 16 band becomes smaller as the excitation energy 
of the relative motion between clusters becomes larger. 
For example, compared to the pure 16 0+ 16 bands, the 
+ states of N=24, 26 and 28 bands are lowered by about 
12, 6 and 2 MeV, respectively. This decrease of the en- 
ergy gain means the enhancement of the 16 0+ 16 molec- 
ular structure in the higher excited bands. Indeed, the 
sum of the w J of the fragmented band members drasti- 
cally increases in N=26 and 28 bands compared to that 
of the N=24 band; they amount to 0.71 and 0.73 for the 
case of the + states of the N=26 and 28 bands. 

In the framework of the deformed-base AMD+GCM, 
the basis states for the GCM calculation are obtained 
from the variational calculation. Since the variational 
calculation optimizes mainly the lowest N=24 band mem- 
bers in which the 16 0+ 16 molecular structure is dis- 
torted, these basis states may be inappropriate to de- 
scribe the N=26 and N=28 bands in which the 16 0+ 16 
molecular structure is drastically enhanced. In other 
words, the present calculation may underestimate the 
enhancement of the 16 0+ 16 molecular structure in 
N=26 and N=28 bands. Therefore, we have included 
the 16 0+ 16 Brink wave functions in the basis states of 
the GCM calculation in addition to the deformed-base 
AMD wave functions obtained from the variational cal- 
culation. The obtained results of the enlarged GCM cal- 
culation are shown in Fig. |3| The enhancement of the 
16 0+ 16 molecular structure has become more promi- 
nent in N=26 and N=28 bands. It is reasonable that 
the excitation energies and the overlaps w J of the N=24 
band members do not change much, since in this band, 
the 16 0+ 16 molecular structure is distorted and the 
inclusion of the pure 16 0+ 16 configuration is less im- 
portant. On the contrary, in the N=26 and 28 bands, the 
excitation energies are lowered by about a few MeV and 
the amounts of the 16 0+ 16 component have increased 
drastically. The sums of the fragmented w J for N=26 and 
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TABLE II: The excitation energies E x and the amounts w J 
of the 16 0+ 16 components of the fragmented + states of 
N=26 and 28 bands. 





N= 

E x 


=26 

w J 


N= 

E x 


=28 

w J 


fragment I 


23.8 


0.54 


31.2 


0.32 


fragment II 


24.0 


0.13 


34.0 


0.45 


fragment III 


25.3 


0.23 


38.7 


0.20 


Eav and ^ w J 


24.2 


0.90 


33.7 


0.98 



28 bands are 0.90 and 0.98, respectively. It is interesting 
that the highest N=28 band members have almost pure 
16 0+ 16 molecular structure when we sum up the frag- 
mented states, and it looks plausible that this band mem- 
bers are assigned to correspond to the observed molecular 
resonances of 16 0+ 16 as is proposed in Ref.[2|. 

Finally, we discuss the fragmentation of the N=26 and 
28 band members. The judgement whether a fragment 
belongs to the N=26 or 28 band is determined by eval- 
uating the amount of the overlaps with the N=26 and 
28 band members obtained from the GCM calculation 
with the pure 16 0+ 16 wave function. As an example, 
the energies and the amounts of the 16 0+ 16 compo- 
nent for + fragments are listed in Table [H] The + 
states of the N=26 and 28 bands are fragmented into 
three states. In the present calculation, the number of 
the fragments do not strongly depend on the angular mo- 
mentum and the principal quantum number N. At most, 
the 12 + state of N=28 band is fragmented into five states. 
These fragmentations are mainly caused by the coupling 
with the states with medium deformation which appear 
as a small peak between normal deformed states and the 
superdeformed states and also by the coupling with the 
16 0+ 16 0* states which are included in the largely de- 
formed (/3 > 0.9) wave functions where 16 0* stands for 
distorted 16 0* cluster. Details of these couplings are im- 
portant to compare the present results and the experi- 
ments and will be investigated in the future. 

To summarize, we have shown that the superdeformed 
band obtained from the HF(B) calculations and the 
Pauli-allowed lowest N=24 band of the 16 0+ 16 molecu- 
lar bands are essentially identical. In this band, 16 0+ 16 
molecular structure is distorted by the effects of the de- 
formed mean-field formation and the spin-orbit force. 
This distortion is not small and lowers the excitation en- 
ergy largely, but this band members still have the consid- 



erable component of the 16 0+ 16 molecular structure. 
We have obtained two excited bands which are generated 
by the excitation of the relative motion between two 16 O 
clusters contained in the 16 0+ 16 component of the su- 
perdeformed band. In the excited N=26 and 28 bands, 
the distortion is less important and the band members 
have the prominent molecular structure. The members 
of these bands are fragmented into several states and the 
assigment of the N=28 band members to the oberved 
16 0+ 16 molecular resonances looks plausible. 
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